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We examine the formation of vortices during tfie nonequilibrium relaxation of a high-temperature 
initial state of an Abelian-Higgs system. We equilibrate the scalar and gauge fields using gauge- 
invariant Langevin equations and relax the system by instantaneously removing thermal fluctuations. 
For couplings near critical, Kc = \/\/e = 1, we observe the formation of large clusters of like-sign 
magnetic vortices. Their appearance has implications for the dynamics of the phase transition, for 
the distribution of topological defects and for late-time phase ordering kinetics. We offer explanations 
for both the observed vortex densities and vortex configurations. 
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Recently much progress has been made in the study 
of the dynamics of phase transitions. The emerging un- 
derstanding may prove relevant for a variety of exper- 
iments, ranging from the collisions of heavy nuclei at 
RHIC PI] t o the sudden cooling of condensed matter sys- 
tems ffl^. Nonequilibrium behavior is particularly im- 
portant for understanding the population of topological 
defects which remains after a transition [||j9j. 

While most previous research focused on global sym- 
metries, it is useful to explore systems with (local) gauge 
symmetry [ |To[ . This includes critical phenomena in the 
early universe, heavy-ion collisions and superconductors 
0. With gauge symmetries, disorder is more subtle as 
local phase gradients can be removed by gauge transfor- 
mations Moreover, gauge fields contribute to the 
(thermo)dynamics of the system, potentially leading to 
new phenomenology. 

In this work we examine how critical dynamics is al- 
tered in the presence of gauge fields. We choose a simple 
and controlled physical setting, describing the behavior of 
long wavelength fields in the quench of 2D superconduc- 
tors. Disorder in the superconducting state is associated 
with the presence of topological configurations (vortices). 
We analyze the formation and dynamics of these defects 
and trace their origins to thermal electromagnetic fluc- 
tuations. The most important result of this work is the 
detection of large clusters of like-sign vortices which form 
in the wake of a thermal quench. 

Hindmarsh and Rajantie (HR) recently argued 
that gauge fields alter the standard Kibble-Zurek predic- 
tions for the distribution of topological defects formed 
in a phase transition. They also observed that magnetic 
fluctuations can create vortex clusters. We provide ev- 
idence strongly correlating the density and distribution 
of defects to thermal magnetic fluctuations. In distinc- 
tion to HR, we argue that an important role is played 
by the critical magnetic fields of a superconductor and 



we exhibit the importance of magnetic fluctuations on a 
variety of large length scales in seeding vortex clusters. 

We investigate the dynamics of a complex scalar field 
minimally coupled to electromagnetism, the Abelian- 
Higgs model, in two dimensions. This provides, for exam- 
ple, a phenomenological description of a superconducting 
film. To accurately compare our results with experimen- 
tal films the effects of a small third dimension need to 
be carefully considered. Nonetheless, this system pro- 
vides a tractable and interesting model, allowing for a 
detailed examination of the importance of magnetic fluc- 
tuations on critical dynamics and defect formation. The 
Lagrangian density is 

C^-\f,^F^'^ + \\D,4>\'-^{\^''-v^)\ (1) 

where cf) is a, complex scalar field, F^i, — i9^A^ — d^A^j^ is 
the field strength tensor for electromagnetic gauge poten- 
tial Af^ and Z?^ = dp, — ieA^,. Of the three parameters in 
the model (A, e and v) two may be rescaled away through 
coordinate and field redefinitions cf) — > (pv, Ap, — s- A^v, 
X vVXx, t v^/Xt and /? = v'^/S, with /3 = (fcsT)-\ 
where ks is Boltzmann's constant. The dimensionless 
parameter k — (the ratio of the London penetra- 
tion depth to the scalar correlation length) controls the 
relative strength of the gauge interactions. With these 
choices the crossover between type-I {n < 1) and type- 
II (k > 1) superconductivity occurs at critical coupling, 
Kc — I. When K — > cx) the global phase invariance of 
the theory is recovered. In the following we work in tem- 
poral gauge, Aq — 0, and therefore obtain Gauss' law, 
V ■ E — 2e Im [0*7r], as a constraint on the evolution. 

To thermalize the system the scalar and gauge fields 
are coupled to a reservoir at a temperature T. This is 
achieved by adding dissipative and stochastic Langevin 
sources to the equations of motion. These cannot break 
gauge invariance, i.e. the Gauss constraint must be pre- 
served, resulting, generally, in multiplicative noise and 
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dissipation terms jlj]. There is more than one gauge in- 
variant form for the equations of motion and we choose to 
couple our system to a reservoir through the observables 
{|</)p, E}. The details of this choice are irrelevant to the 
state of canonical equilibrium reached by the system at 
long times. 

After equilibrating the system, we instantaneously re- 
move the noise and evolve with overdamped dynamics 



dt(j>a = 



1 



Ji = -e^l^pAj - eeab4'adi(f>b. 



(2) 



The system of equatio ns (p|) , expressed in terms of lattice 
gauge invariant fields |l5| , was solved numerically. The 
Gauss constraint was satisfied to about 0.1%. We used 
lattice spacing dx — 0.5 and time step dt — 0.02. 

A snapshot taken after the quench reveals striking 
structure in the spatial distribution of vortices. In Fig. |l| 
we provide a contour plot of the magnetic flux across 
a section of the lattice with e = 0.5 after a quench from 
T — 10. The critical temperature is Tc ~ 0.35. The snap- 
shot is late enough in the quench that the magnetic flux 
is almost entirely localized into vortices, each of which 
carries the fundamental flux quantum $o = The 
large clusters of like-sign vortices are an obvious depar- 
ture from the spatial vortex distribution observed after 
a quench of purely scalar fields. In the scalar case, de- 
fects and antidefects appear primarily as widely sepa- 
rated and anticorrelated pairs [|l8|. The appearance of 
closely-packed vortices is also surprising because (though 
they are not strongly coupled) like-sign vortices repel. 

To explain the clustering of vortices we explore the 
dynamics of the electromagnetic fields. In equilibrium, 
at sufficiently high T above the transition, {BkB-k) — T. 
When we remove the external noise, all fields dissipate. 
However, not all length scales dissipate at the same rate. 
Long wavelength magnetic fluctuations suffer an analog 
of critical slowing down due to the absence of a magnetic 
thermal mass in an Abelian gauge theory |p^ . Under 
these circumstances we expect 



Bk{t) ^ Bk{G)e^j>[-kH] . 



(3) 



While strictly true only for free dynamics, we verified nu- 
merically that Eq. (^) adequately describes the behavior 
of B in the coupled system for the initial stages of cooling. 
This overdamped evolution results in the freezing of mag- 
netic fluctuations on large spatial scales. Fluctuations in 
the scalar field are also damped during the cooling. When 
fluctuations of the scalar field have decayed sufficiently, 
long-wavelength modes feel the spinodal instability and 
the scalar field rolls to its nonzero minimum. 




FIG. 1. Contour plot of the magnetic flux across a sec- 
tion of tiie lattice after the quench. White, localized, regions 
denote vortices while black regions denote antivortices. The 
giant like-sign vortex clusters are surprising and are a result of 
the quench dynamics in a type-II superconductor, beginning 
with the formation of superconducting regions (in which the 
gauge fields are small B < Bc\) surrounding normal islands 
of coherent field. As the quench proceeds the normal islands 
fragment into clusters of individual vortices. 



The evolution of the system can be qualitatively under- 
stood with reference to the state of minimal free energy in 
an external magnetic fleld. There are three possibilities 
|T7f delimited by 



Bel = $0 log K/47reG, B,2 = *o/2<s. 



(4) 



where ^5 and S^q are the scalar and gauge correlation 
lengths. Regions where B > Bc2 are forced to remain in 
the normal phase. For Bd < B < Bc2 the system enters 
a vortex state, where B penetrates a superconducting 
state in an array of like-sign vortices. For B < Bd the 
magnetic field is completely expelled and the system be- 
comes a spatially homogeneous superconductor. In the 
presence of large thermal B, the quench dynamics cause 
the system to enter the vortex state in certain regions of 
space. Once formed the vortex state is long-lived because 
topological charge is conserved and the evolution is over- 
damped. Where B < Bd, no vortices are formed by this 
mechanism; the flux is expelled by a Meissner effect and 
the picture of defect formation in terms of the dynamics 
of the scalar fleld alone is recovered. We arrive at a pic- 
ture of spontaneous vortex formation with two primary 
qualiflcations: 

1. The number of vortices formed in the transition is 
proportional to the amount of magnetic flux present 
when the system develops supercurrents. 

2. The magnetic field splinters into vortex clusters 
when B > Bd- The clustered vortex distribu- 
tion reflects the long-range correlations of the ini- 
tial magnetic fleld, amplifled by the expulsion of 
flux from regions where B < Bd- 
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Prediction 1 was first made by HR and our analysis of- 
fers quantitative evidence that this is correct. However, 
our second prediction is based on a quahtatively differ- 
ent picture of the dynamics of the gauge fields and order 
parameter during the phase transition. 

To test the first prediction we examine the magnetic 



flux, |ci>| = / 
value then 



'B^(Px. If we replace by its average 



|$| - V<B^>A, 



(5) 



where A is the total area of the film. The determination 
of the time when supercurrents develop and therefore the 
point at which to evaluate \J <B'^> depends on the cou- 
pled dynamics of the scalar and electromagnetic field and 
will be addressed elsewhere. However, we can make the 
following (parametric) prediction ||l3|] , 



Pa. 



,f = Ce\/T, 



(6) 



where pdcf is the (area) density of vortices and T is 
the initial temperature. C carries dimensions of inverse 
length and parameterizes the details of the cooling dy- 
namics. In practice we found C 0.01 demonstrating 
that only a small fraction of the original thermal mag- 
netic energy is available for vortex formation. We expect 
Eq. (^ to hold only when the magnetic field at the time 
of formation is greater than Bci ■ 




FIG. 2. a) pdef X 10 vs. e from numerical evolutions. The 
solid line is the best linear fit. b) pdcf x 10^ vs. T from 
numerical simulations. The solid line is the best power-law 
fit p ~ T", where a = 0.43 ± 0.02. In both plots error bars 
denote standard deviations among stochastic realizations. 

Figs. I^a and |^b show pdef vs. e and the initial tempera- 
ture T respectively. Both are in agreement with Eq. (|^). 
The slightly smaller exponent in Fig. ||b is an indication 
that some fluctuations decay before vortices are formed. 




FIG. 3. a) Nc{l) for the cluster distribution of Fig. |l| (cir- 
cles) compared against a random distribution with the same 
number of defects (triangles). The solid line shows the be- 
havior expected from equilibrium fluctuations in the mag- 
netic flux, "l>net = VtL The close agreement between Nc{l) 
measured from the vortex distribution and Nc{l) calculated 
from equilibrium is strong evidence that vortex clusters in the 
gauge quench originate in equilibrium magnetic fluctuations, 
b) Nc{l) for the purely scalar theory, e = 0, normalized to the 
same density as Fig. 3a. The solid line is the best fit to \/I 

While the number of vortices depends on the details 
of the cooling dynamics, the clustered distribution of 
vortices reflects initial magnetic thermal fluctuations on 
many scales. During the quench, magnetic fields are first 
damped and then redistributed by the Meissner effect: 
expelled where B < Bd into initially normal magne- 
tized islands, which eventually fragment into vortices. 
To quantify the properties of the vortex cluster distri- 
bution we examine the average excess topological charge 
Nc{l) within a square of side /. Nc{l) can be measured 
in experiments using e.g. a SQUID, sensitive to flux over 
an area P. To calculate Nc{l) we compute the absolute 
value of the net topological charge in a box of side length 
I centered on a lattice point. We then average over ev- 
ery lattice point. Fig. ^ shows the result for a typical 
run with e = 0.5. Nc{l) observed for the gauge quench 
is indicative of clustering and is larger than that of a 
random distribution of the same number of defects. At 
small / the presence of clusters of like-sign vortices leads 
to faster than linear growth of Nc{l). At large I both 
distributions increase linearly. This is expected when de- 
fects are positioned at random: with average density pdof 
fluctuations are of order y/ N{1) where N{1) is the average 
number of defects in box of size /. Therefore 



N- 



random 



(0 



(7) 



The linear behavior and larger slope of Nc{l) measured 
from a gauge quench suggest that clusters of like-sign 
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defects are randomly distributed. 

In scalar systems defect distributions are not random 
psf . Fluctuations of the net topological number over an 
area occur due to the fluctuations of the field phase along 
the boundary. After a quench, the phase is correlated on 
a scale ^ with 1/$^ separate domains along a path of length 
/. Thus, 



NS) - (Pdcf)'/'VI. 



(8) 



Fig. ||b shows Nc{l) measured from simulations of a 
purely scalar theory. Eq. (^) has also been verified exper- 
imentally in nematic liquid crystals ||l9[] and is roughly 
consistent with measured fluctuations in the net flux 
trapped within a loop of Josephson junctions §. 

It is instructive to examine defect formation in the 



presence of an imposed magnetic field (see also 20 1). We 
prepare an uncoupled initial state in which the scalar 
field is thermalized to a temperature close to Tc and the 
magnetic field is 



B{t = 0) = Bo cos ( — 



(9) 



At the time of the quench we turn on the coupling (to 
e = 0.5) and evolve with the overdamped equations (y). 
The results for two different Bq are shown in Fig. m 




FIG. 4. Contour plot of the magnetic flux after a quench in 
an external magnetic field. The left plot corresponds to am- 
phtude Bo ~ 0.01. The right plot is for amplitude Bo ~ 0.16. 
The lower critical field is Bci = 0.08 while Bc2 = 0.5. 



For small amplitudes, the final defect distribution is not 
heavily influenced by the magnetic field and large vortex 
clusters are absent. In contrast for B > Bd the magnetic 
field has a profound influence on the final defect distri- 
bution. These simple considerations mimic the behavior 
occurring on many length scales in the thermal quench. 

In conclusion, we have shown that the critical dynam- 
ics of 2D superconductors differs substantially from scalar 
field models. Specifically we observed that for quenches 
from high T and for remnants of long-wavelength 

magnetic fields seed the formation of topological defects 
in much larger numbers and in strikingly different con- 
figurations than scalar field theories. We established the 
dependence of the density of vortices on e and the initial 



temperature T and we connected the large-scale prop- 
erties of the vortex distribution to the spectrum of ini- 
tial thermal magnetic fluctuations. We also argued for 
the importance of Bd in determining the effect of mag- 
netic fluctuations upon defect formation. With recent ex- 
perimental progress in vortex imaging pl[ | , clusters such 
as Fig. |l] might be discernible in superconductors with 
K ~ Kc (e. g. Nb). A detailed analysis of the statistical 
properties of these clusters is under investigation. 

We thank M. Hindmarsh, A. Rajantie and N. Golden- 
feld for useful comments. This work was partly supported 
by D.O.E grant DF-FC02-94ER40818. Numerical work 
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